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Using functional renormalization group method, we study the favorable condition for electronic 
correlation driven time reversal invariant topological superconductivity in symmetry class DHL For 
non-centrosymmetric systems we argue that the proximity to ferromagnetic (or small wavevector 
magnetic) instability can be used as a guideline for the search of this type of superconductivity. This 
is analogous to the appearance of singlet unconventional superconductivity in the neighborhood of 
antiferromagnetic instability. We show three concrete examples where ferromagnetic-like fluctuation 
leads to topological pairing 
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Topological insulators and superconductors have be- 
come a focus of interest in condensed matter physics. [1, 2] 
These states are characterized by symmetry protected 
gapless boundary modes. The existence of these modes 
reflects the fact that it is impossible to deform a topo- 
logical insulator/superconductor into its non-topological 
counterpart without crossing a quantum phase tran- 
sition. The free-fcrmion topological superconductors 
and insulators have been classified into ten symmetry 
classes [3, 4]. In each space dimension precisely five of 
these classes have representatives. Examples of topo- 
logical insulators include the T-brcaking integer quan- 
tum Hall insulator (2DEG[5]), and the T-preserving 
topological insulators in two and three space dimen- 
sions (2D HgTc quantum wells [6] and 3D Z 2 topolog- 
ical insulators [7]). Examples of topological superfluid 
or superconductor include the T-breaking 3 He-A[8] and 
Sr 2 Ru0 4 [9] (likely), and T-invariant 3 He-B[8]). 

In this fast growing field discovering new topological 
materials is clearly one of the most important tasks. 
In this regard predicting topological superconductors is 
much harder than predicting topological insulators. This 
is because knowing the desired Bogoliubov de Gennes 
(BdG) band structure[10] only meets half of the chal- 
lenge. The other half requires the knowledge the mi- 
croscopic interactions which favor the desired quasipar- 
ticle band structure as the mean-field theory. There 
arc many interesting proposals for inducing topologi- 
cally non-trivial superconducting pairing via the prox- 
imity effect. [11-13] In these proposals, pairing is arti- 
ficially induced by a (non-topological) superconductor. 
The reason the induced superconducting state is topo- 
logical is due to the novel spin-orbit coupled electronic 
wavefunctions in the normal state. A notable exception 
is the intriguing proposal that the superconducting state 
of Cu^Bi2Se3 is topological. [14] 

Leaving topology aside, it is extremely challenging to 
predict superconductivity itself. This is because the en- 
ergy scale involved in Cooper pairing is usually much 



smaller than the characteristic energies in the normal 
state. However in the last five years various types of 
renormalization group methods have been used to com- 
pute the effective interaction responsible for the Cooper 
pairing in iron-based superconductors. [15] They lead to 
a proposal for why the pairing scale of the pnictides is 
high: the scattering channels triggering antiferromag- 
netism and Cooper pairing have overlaps. Through these 
overlaps strong antiferromagnetic fluctuation enhances 
superconducting pairing. [16] This is consistent with the 
widely known empirical fact: strong superconducting 
pairing often occurs when static antiferromagnetism dis- 
appears. 

There are two classes (Dili and CI) of time-reversal- 
invariant (TRI) topological superconductor in three 
dimensions. [3, 4] They are differentiated by the trans- 
formation properties with respect to time reversal and 
particle-hole conjugation. In this paper we will focus 
on the so-called class Dili, for it has realization in all 
space dimensions. We ask "under what condition is 
time-reversal symmetric topological superconductivity fa- 
vored?" We shall argue that it is near the ferromagnetic 
(to be precise small wavevector magnetic) instability. 
However due to the exponential growth of computation 
difficulty with space dimensionality we shall limit our- 
selves to two dimensions. We shall also restrict the dis- 
cussion to systems with a special type of spin-orbit inter- 
action - the Rashba coupling. This type of spin-orbit cou- 
pling breaks the parity symmetry, hence the supercon- 
ductors under consideration are non-centrosymmetric. 
For discussions of topological pairing in centrosymmet- 
ric systems see, e.g., Rcf. [17]. Many real superconduct- 
ing materials are non-centrosymmetric. Examples in- 
clude CePt 3 Si[18], CeRhSi 3 [19], CeIrSi 3 [20], and the su- 
perconductivity found at the interface of LaA10 3 and 
SrTi0 3 [21]. 

Here are our main results. We present three differ- 
ent mechanisms for topological pairing. We warm up by 
studying a one band model mimicking strongly correlated 
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FIG. 1: A generic 4-point vertex Ti234 is rearranged into P-, 
C-, and D-channels in (a)-(c), respectively. Here k, q, p are 
momenta, /i, v, a, X are spin indices, and m, n denote the form 
factors. On each side of the diagrams, the spin (and sublat- 
tice) labels are absorbed into the form factor labels wherever 
applicable (see the main text). 



fermions in continuum. As we shall discuss this model is 
relevant to the topological superfluidity in the B phase of 
3 He. Next we discuss two other different mechanisms for 
topological pairing. In each case there is a finite param- 
eter range where triplet pairing occurs in the presence 
of ferromagnetic (or small wavevector magnetic) fluctu- 
ation. We explain why, under such condition, a small 
Rashba coupling can induce topological superconductiv- 
ity. We also explain why topological pairing does not 
happen in singlet-dominated materials. The paper con- 
cludes with a guideline for the search of TRI topological 
superconductivity in non-centrosymmetric systems. 

Method - Technically this work requires us to gener- 
alize the functional renormalization group (FRG) ap- 
proach [15, 16, 22, 23] to Hamiltonians without spin 
rotation symmetry. In addition, because the necessity 
to study small momentum transfer particle-hole scat- 
terings we use a Matsubara frequency rather than mo- 
mentum cutoff. All calculations are carried out using 
the singular-mode functional renormalization group (SM- 
FRG) method. [23, 24] 

Consider a generic fully-antisymmctized irreducible 
4-point vertex function Ti234 in 'i'l^i^ ri234) , f / 3 , I'4. 
Here 1,2,3,4 represent momentum and spin (and 
sub-lattice) indices. Figs.l(a)-(c) are rearrangements of 
Ti234 into the pairing (P), crossing (C) and direct (D) 
channels each characterized by a collective momentum 
q. In each channel the vertex function is decomposed 
as Eq. (4) of the Appendix. There {/,„} is a set of 
orthonormal lattice form factors. [25] The spin (and 
sublatticc) indices are contained in the label of the form 
factors as shown in Figs.l(a)-(c). The decomposition in 
Eq. (4) is exact if the form factors are complete, but 
a few of them are often enough to capture the leading 
instabilities. [23, 24] The FRG flow equations for P,C 
and J) as a function of the cutoff scale A are given 
by Eqs.(6),(7) and (8) of the Appendix. The effective 
interaction in the particle-particle (pp) and particle-hole 
(ph) channels are given, respectively, by V pp = —P and 
Vph = C. [Because of antisymmetry D (= — C) does 
not yield any new information.] During the FRG flow 
we monitor the singular values of the matrix functions 



V pp / p h{c\). The most negative singular values, S pp / p h, 
occur at special momenta q pp / p h- While q pp is usually 
zero, <\ p h can evolve under RG before settling down to 
fixed values. The eigen function associated with S pp 
is used to construct the gap function. More technical 
details can be found in the Appendix. 

A strongly correlated one-band model in con- 
tinuum limit - We consider spin- 1/2 fermions hopping 
on a square lattice. The Hamiltonian is given by 

H = ]T*jje(k)a + A7(k) 

k 

+ u Y n n n n + V Y1 n * n J = h o + Hl (l) 

Here ^ = (ipp e(k) = — 2t(cos k x + cosk y ) — [i (t is 
the nearest neighbor hopping integral and fi is the chemi- 
cal potential), i labels the lattice sites, riia = tp^tpia and 
m = rii a . In addition, <jq is the 2x2 identity matrix 
and a denotes the three Pauli matrices. In the following 
we shall set U = 8t and V = — 2t for the on-site and 
nearest neighbor interactions. For the Rashba spin-orbit 
coupling we consider 7(k) = (— sin k y , sin k x , 0). 

Combining the time-reversal and point group (C^y in 
the present case) symmetries , it can be shown that 
the Cooper pair operator = ^ k 'f^Ak* 1 ^ takes the 

form, [26] A(k) = [0(k)ero + d(k)-cr]io-2, where d(k) trans- 
forms, under the point group, like the product of 0(k) 
and 7*(k). In the cases we have studied, to a good ap- 
proximation, we can write 

A(k) = [^(k)ao+x(k) 7 (k)-o ! ]i C 7 2 , (2) 

where j(k) = 7*(k)/|7(k)|, </>(k) and %(k) are even func- 
tions of k (real up to a global phase) and transform 
according to the same irreducible representation of the 
point group (for multi-dimensional representations there 
are several <f> and x' s )- I n Landau theory, </> and x a °t 
as order parameters, and can induce each other in the 
presence of the Rashba coupling (A ^ 0). 

It is important to note that the Rashba term splits 
each of the otherwise spin-degencrate Fermi surface into 
two. The spin split Fermi surfaces are characterized by 
eigen values ±1 of j(k) ■ a. In the case where 4>(k) and 
%(k) are nodeless, the gap function on the two split Fermi 
surfaces will have opposite sign if the magnitudes of x(k) 
dominates over </>(k). It turns out that for each pair of 
Fermi pockets surrounding a TRI k point the above sign 
reversal leads to two counter-propagating Majorana edge 
modes . Thus topological pairing requires the triplet \~ 
component to be dominant. Moreover sign reversal (in 
the gap function) on an odd/even pairs of the spin-split 
Fermi surfaces (satisfying the condition specified above) 
will lead to strong/weak topological superconductivity. 
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For A = O.Olt and fi = —M the spin-split Fermi sur- 
faces are shown in Fig. 2(a). The pockets are small, mim- 
icking the continuum limit. The form factors used in our 
SM-FRG extend up to second neighbors in real space. [25] 
The RG flow of Spp/p^ are shown in Fig. 2(b). During 
the flow (\ p h evolves from qi = (n, n) and settles down 
at q2 = 0. By inspecting the spin structure of the q2- 
singular mode we find it corresponds to ferromagnetic 
fluctuation. The leading pairing channel is extended s- 
wave at cutoff energies above point A (because the bare 
V < 0). But at lower cutoff energies the increased ferro- 
magnetic fluctuation around q2 enhances pairing in the 
triplet channel via their mutual overlaps (sec Appendix). 
The cusps at A, B and C associated with the S pp flow is 
due to the evolution of the leading pairing form factor. 
The gap function is determined by the singular mode as- 
sociated with S pp at the diverging cutoff scale. The result 
is a dominant x-component together with a much smaller 
</>-component. The corresponding gap function on the 
two Fermi surfaces is shown in Fig. 2(a) (gray scale). A 
sign change is clearly visible. According to the estab- 
lished criterion, [2] this pairing state is topological. To 
verify this, we calculate the BdG energy spectrum us- 
ing the obtained pairing form factor in a strip geometry 
(open-boundary along x). The resulting eigen energies 
as a function of q = k y is shown in Fig. 2(c). There are 
two in-gap branches of Majorana edge modes associated 
with each edge. 

Had we turned off the Rashba coupling, the leading 
pairing channel (p-wave) would be two-fold degenerate 
(with dominant amplitudes on 1st neighbor bonds). 
Under this condition even an infinitesimal Rashba 
coupling breaks the degeneracy by linearly recombin- 
ing the p- waves into A(k) = i sin k x ao + shikya^, or 
x(k) = |-y(k) I in Eq.(2), leading to a gap function ±x(k) 
on the infinitesimally split Fermi surfaces. This gap 
function has the same symmetry as the two dimensional 
version of the 3 He B phase. In fact there are strong 
similarities between the model (and its properties) 
described above and the B phase of 3 He. For example 
The small filling fraction enables this model to describe 
the continuum limit. The strong on-site repulsion 
mimics the short-range strong repulsive correlation in 
3 Hc, while the weaker nearest neighbor attraction is 
a caricature of the tail of the Lennard- Jones poten- 
tial. The enhanced ferromagnetic fluctuation and the 
resulting triplet pairing is consistent with the pairing 
mechanism described by Anderson and Brinkman[27]. 
In addition, the Rashba coupling plays a similar 
role as the parity -invariant spin-orbit interaction in 
3 He: they both lift the degeneracy in the pairing channel. 
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FIG. 2: (Color online) Left panels: Results for the toy model 
of He 3 , (a) The Fermi pockets and the associated gap func- 
tions (gray scale, in units oii). The spin-splitting between the 
pockets is enlarged for clarity. The box is the zone boundary, 
(b) The SM-FRG flow of S pp / P h versus cutoff scale A. Arrows 
mark sharp changes in the RG evolution, (c) The low energy 
BdG eigen spectrum in a strip (open along x) as a function of 
the momentum qy. Right panels: the same plots as in the left 
panels but for a model whose Fermi surfaces are in proximity 
to van Hove singularities. 



Here we set V = (so that the interaction is purely 
repulsive), and add a 2nd neighbor hopping t' so that 
e(k) = — 2t(cosfc x + cos k y ) — At' cos k x cos k y — /i. For 
t' = -0.475t, \i = -2t and A = 0.0H, the spin-split 
Fermi surfaces are shown in Fig. 2(d). They are pointy 
along x and y reflecting the existence of saddle points 
(van Hove singularities) on the Brillouin zone boundary. 
These features lead to enhanced ferromagnetic corre- 
lations via the Stoner mechanism. As a consequence 
triplet pairing is enhanced and the diverging scale of 
Spp is raised. These are shown in Fig. 2(e) for U = 2.5t. 
The arrows associated with the S p h flow record the q p /j 
evolution from qi = (tt, it) to q2 = 0. As in the previous 
section, this implies ferromagnetic correlations at low 
energies. The final gap function on the Fermi surfaces 
are shown in Fig. 2(d) (gray scale), again with desired 
sign reversal. Such a pairing function leads to the BdG 
energy spectrum shown in Fig. 2(f) in a strip geometry. 
The edge modes are apparent. 



Topological pairing in the vicinity of van 

Hove singularity— In this section we demonstrate 

that topological pairing can be enhanced when the 

Fermi surfaces are close to van Hove singularities. 



Topological pairing enhanced by inter-pocket 
scattering — In this section we show a third route to 
topological pairing. In this case pairing is triggered by 
inter- Fermi surface scattering in a way similar to the pair- 
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ing in the pnictides[15, 16]. 

Consider a honeycomb lattice. The single particle 
Hamiltonian is given by 

H Q = ^Its^i+S ~ iA ^ $l(z x S nn ■ a)^ i+Snn 

-HY,*\%- (3) 

i 

Here i labels lattice sites, 8 runs over the 1st and 2nd 
neighbor bonds, with t& = t,t' . The spin-dependent hop- 
ping, the Rashba term, is limited to the nearest neigh- 
bor bonds 8 nn . Choosing a lattice site as the origin, 
the point group is C$ v . For the SM-FRG calculation, 
we choose the form factors up to the 2nd neighbors. [25] 
(Since the honeycomb lattice has two sites per unit cell 
the labels of the form factors in Fig.l include the sub- 
lattice indices. [23].) The Fermi surfaces for t' = 0.357t, 
A = 0.02i and \i = 1.664i are shown in Fig. 3(a). There 
are a few interesting features of the band structure that 
are worth noting (1) The Fermi surfaces encircle either 
the zone center (r) or the zone corners (K and K'). How- 
ever only r is TRI, hence according to Ref. [2] only the 
T-Fermi surfaces are topologically relevant. (2) The T 
and K-pockcts have close by segments, hence allow small 
momentum transfer particle-hole scattering. If such scat- 
tering is magnetic, it corresponds to nearly ferromag- 
netic fluctuations, hence can induce triplet and topolog- 
ical pairing. 

In the following we show for U = 1.26i this is ex- 
actly what happens. During the RG flow shown in 
Fig. 3(b), the strength of S p h increases and q p h evolves 
fromqi = (0.667, 1.152)tt to q 2 = 0, q 3 = (0.250, 0.048)tt 
and finally settles down at q4 = (0.333, 0.192)7r. We 
have checked that q4 corresponds to the scattering be- 
tween near by parallel segments between the T and K 
pockets. Inspection of the spin structure of the singular 
mode associated with q2,3.4 reveals that they corresponds 
to spin fluctuations. As such fluctuations are enhanced, 
they causes S pp to grow in magnitude and eventually di- 
verge at a relatively high critical scale. The resulting 
gap function is shown in Fig. 3(a) in gray scale. It is fully 
gapped on all Fermi surfaces, and have opposite sign on 
each pair of spin-split pockets. Since the K-pockets are 
topologically irrelevant, the sign change between the T 
Fermi surfaces implies the pairing is strong-topological. 
To verify this we consider a strip schematically shown in 
Fig. 3(c). It is open along a and periodic along b direc- 
tions. The BdG energy spectrum as a function of the 
momentum q = k ■ b is shown in Fig. 3(d). There are two 
branches of Majorana modes at each edge. 

Thus in all of the above examples we have seen small 
momentum magnetic fluctuations degenerate triplet 
pairing, and degenerate triplet pairing + Rashba cou- 
pling topological pairing. The fact that ferromag- 
netic fluctuations enhance triplet pairing has a long his- 
tory. These include the works on the pairing of 3 Hc,[8], 
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FIG. 3: (Color online) (a) The Fermi pockets and the as- 
sociated gap functions (gray scale, in units of t). The spin- 
splitting between each pair of Fermi pockets is enlarged for 
clarity. The hexagon is the zone boundary, (b) The SM-FRG 
flow of Spp/ P h versus the cutoff scale A. Arrows mark sharp 
changes of q p h during the RG flow, (c) A strip (marked by 
the thick lines) open along a and periodic along b directions 
(a and b are primitive lattice vectors), (d) The low energy 
BdG eigen spectrum for (c) as a function of the conserved 
momentum q along b. 

and the extension of the Kohn-Luttinger theorem to p- 
wave pairing for 2D and 3D electron gas in the dilute 
limit. [28, 29] For lattice systems, a 2D Hubbard model 
with large enough 2nd neighbor hopping has been shown 
to exhibit p-wave pairing for small band fillings. [30] 

Finally, if pairing is singlet in the absence of spin-orbit 
interaction a weak Rashba coupling only induces a small 
triplet component, hence is insufficient to induce the de- 
sired sign change in the gap function. Of course this docs 
not rule out the possibility of topological pairing in the 
presence of strong spin-orbit interaction. 

In conclusion through the study of the above three, 
and many other not included, examples we conclude that 
TRI topological superconductivity in symmetry class 
Dili should occur in systems close to the ferromagnetic 
(or small wavevector magnetic) instability. Bandstruc- 
ture wise, in the absence of Rashba coupling, these 
systems should have an odd number of spin-degenerate 
Fermi pockets (each enclosing a TRI momentum) in 
order for strong topological pairing to occur. 
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APPENDIX 

In the supplementary material, we provide the techni- 
cal details of the SM-FRG method. [23]. 

We begin by reviewing the definition of the vertex func- 
tions used in the main text. Consider a generic fully- 
antisymmetized irreducible 4-point vertex function 1^1234 
in 

1234)^3^4- Here 1,2,3,4 represent momen- 
tum and spin (and sublattice) indices. Figs.4(a)-(c) are 
rearrangements of Ti234 into the pairing (P), crossing (C) 
and direct (D) channels each characterized by a collec- 
tive momentum q. The rest momentum dependence of 
the vertex function can be decomposed as, 

^k+q,-k,-p,p+q — /m(k)-Pmn(q)/rt(p), 

ran 

r k+q A p,k,p+q ~* fm s ) C m n (q) fn (p) ) 

ran 

rCr q A p, P+q ,k /m(k)An«(q)/«(p). (4) 

ran 

Here {f m } is a set of orthonormal lattice form factors. 
The spin (and sublattice) indices are contained in the 
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FIG. 4: A generic 4-point vertex Ti234 is rearranged into P-, 
C-, and D-channels in (a)-(c), respectively. Here k, q, p are 
momenta, /i, v, a, A denote spins, and m, n denote the form 
factors. On each side of the diagrams, the spin (and sublat- 
tice) labels are absorbed into the form factor labels wherever 
applicable. The one-loop diagrams that contribute to dP, dC 
and dD are shown in (d)-(f), respectively. 



label of the form factors as shown in Figs.4(a)-(c). The 
decomposition in Eq. (4) is exact if the form factors are 
complete, but in practice a few of them are often enough 
to capture the leading instabilities. [23, 24] Because of full 
antisymmetry, the matrices C and D satisfy D = — C, 
and are therefore not independent. In the following D is 
used for bookkeeping purpose. 

Ignoring the spin and sublattice labels for the moment, 
the form factors are given by 



/m(k) = /™( r ) exp(-ik ■ r), 



(5) 



where / m (r) transforms according to an irreducible rep- 
resentation of the point group, and r is the bond vec- 
tors connecting the two 'S's (or two vpt's) in Fig. 4(a) 
and one iff and one vpt in Fig. 4(b) and (c). In our cal- 
culation we choose form factors up to the 2nd neighbor 
bonds. We have checked that longer range form factors 
does not change the results qualitatively. To be specific, 
for square lattice, the real-space form factors we used are 
1) f x = 1 for on-site; 2) f 2 = 1/2, f 3 = (l/2)cos20 r , 
f± = -y/l/2cos 9 r , and f$ = y/f/2 sin 8 r for 1st neigh- 



bors, where 6 T is the azimuthal angle of r; 3) = 1/2, 
f 7 = (l/2)sin20 r , /„ = yi72cos(# r - tt/4) and f 9 = 
y/1/2 sin(# r — 7r/4) for 2nd neighbors. For hexagonal lat- 
tices, the form factors we used are 1) fi = 1 for on-site; 
2) f 2 = yi73, h = y/2/3cos6 r and / 4 = y/2/lsme t 
for 1st neighbors; 3) fs — \/l/6, fe — ^/l/3cos 9 T , 
h = y/T/Emn0 r , fs = V /T73cos20 r , f 9 = ^T/3sm26 r , 
fw = \/l/6cosSd r for 2nd neighbors. Notice that the 
1st neighbor bonds stem from different sublattices are 
negative to each other. 

In the case where sublattices are involved, the form fac- 
tor label m also includes the sublattice indices associated 
with the two "f's (or ^''s), or the iff and . However, 
once r is fixed only one of these sublattice indices is in- 
dependent. We include the independent sublattice index 
in the form factor label, (to, a) — > m. Here a labels, e.g., 
the fermion field 1 or 4 in Fig. 4(a), 1 or 4 in (b), and 1 
or 3 in (c). The sublattice index is an independent label 
because point group operations do not mix sublattices 
when the origin is chosen to be a lattice site. 

The total number of form factors TV in a calculation 
is determined by the number of real space neighbors, 
the number of sublattices and the four spin combinations 
(ji, v) =|t, t4-> It, ii associated with two * (P channel) 
or the * and (C and D channels). Thus P, C and D 
are all ./V x N matrix functions of momentum q. 

The Feynman diagrams associated with one-loop con- 
tributions to the flow of the irreducible 4-point vertex 
function are given in Fig.4(d)-(f). They represent the 
partial changes dP, dC and dD, respectively. (Notice 
that the three diagrams in Fig.4(d)-(f) become the usual 
five diagrams in the spin-conserved case.) The internal 
Greens functions are convoluted with the form factors 
hence in matrix form, 



dP/OA = P X ' Pp P/2, 

oc/dA = c x ' ph c, 

dD/dA = -D x ' ph D, 



(6) 



where we have suppressed the dependence of the collec- 
tive momentum q, and 



(Xp P ) mn = ^ J ^ J ^^f m (p)G(P + <l,^n)G(-p,-iuj n )f*(p)e(\uj n \ - A) 
= ~ [ P^-fm(p)G(p + q, *A)G(- P> -iA)/*(p) + (A -> -A), 

27T J bBZ 

(X' ph ) mn ^-^j^j J~"/m(P)C(P + «3) iu} n)G(p, iuj n ) f*(p)6(\uj n \ - A) 

= -5- / ^-/m( P )G(p + q, *A)G(p, iA)/:(p) + (A -> -A), (7) 
2tt J S BZ 
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where G is the free fermion Greens function, Sbz is the 
total area of the Brillouine zone. Here A > is the in- 
frared cutoff of the Matsubara frequency uj n . As in usual 
FRG implementation, the self energy correction and fre- 
quency dependence of the vertex function are ignored. 

Since dP, dC and dD come from independent one- 
loop diagrams, they contribute independently to the full 
<£Ti234, which needs to be projected onto the three chan- 
nels. Therefore the full flow equations are given by, for- 
mally, 

dP/dk = dP/dk + P(dC/dk + dD/dk), 
dC/dk = dC/dk + C(dP/dk + dD/dk), 
dD/dk = dD/dk + D(dP/dk + dC/dk), (8) 

where P, C and D are the projection operators in the 
sense of Eq. (4). Here we have used the fact that 
K(dK) = dK for K = P, C, D. In Eq. (8) the terms 
preceded by the projection operators represent the over- 
laps of different channels. For two channels to overlap, 
the spatial coordinates of all four fermion fields must lie 
within the range set of the form factors. In the actual cal- 
culation the projections in Eq.(8) are preformed in real 
space. 

The effective interaction in the particle-particle (pp) 
and particle-hole (ph) channels are given, respectively, 
by Vp P = — P and V p h = C. By singular value decompo- 
sition, we determine the leading instability in each chan- 
nel, 

Vr(lx) = J2 S % ( t>x(rnm(n), (9) 

a 

where X = pp,ph, S x is the singular value of the a-th 
singular mode, (j) x and ip x are the right and left eigen 
vectors of Vx , respectively. We fix the phase of the eigen 
vectors by requiring Re[J^ (f> x (m)ip x (m)] > so that 
S x < corresponds to an attractive mode in the X- 
channcl. 

In the pp-channel q pp = corresponds to the zero 
ccnter-of-mass momentum Cooper instability. The ma- 
trix gap function A k in the spin and sublattice basis is 
determined as follows. A singular mode </>" leads to a 



pair operator (in the momentum space), 

* k Ak*l T k = E ^(k)^(m)/ m (k)>t m(y( _ k) . 

(10) 

Here a is the independent sublattice index, and a m is 
the second sublattice index determined by a and m as 
discussed earlier, and [i, v are spin indices. The parity 
of At under space inversion determines the singlet and 
triplet components. The gap function in the band eigen 
basis can be determined by the unitary transformation 

% = *M, (ii) 

where the columns of are the Bloch states {|k, n)} (n 
is the band index). Under Eq. (11) the pairing matrix 
transforms into 

A k = (7 k A k [/^ k . (12) 

In the weak coupling case (i.e., when the magnitude of 
the superconducting gap is much smaller than the band- 
width), only the diagonal part of A (i.e., intra- Fermi sur- 
face pairing) is important. Since Eq. (12) involves Bloch 
states at two different momenta, the phases of the asso- 
ciated Bloch states enters A. Since there is time- reversal 
symmetry we fix the Bloch state phase at k and — k by 
demanding T|k, n) = | — k, n) and T 2 |k, n) = — |k, n), 
where T = ia%K is the time-reversal operator. 

In the particle-hole channel, we calculate the singular 
values associated with V p h(<\) at all momenta q. Unlike 
the Col[oooper channel, the most negative singular value 
can occur at non-zero momentum q p /j. The associated 
particle-hole operator is given by 



E i(k+<i)ftwr m (k^,(k). (is) 

m—(m,a,fi,i/) 

Usually the on-site form factor dominates in the particle- 
hole channel. By inspecting the spin structure of the 
on-site form factor one can easily determine whether the 
instability is charge or spin like. 



